Nonlinear optical probe of tunable surface electrons on a topological insulator 
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We use ultrafast laser pulses to experimentally demonstrate that the second-order optical response 
of bulk single crystals of the topological insulator Bi2Se3 is sensitive to its surface electrons. By 
performing surface doping dependence measurements as a function of photon polarization and sam- 
ple orientation we show that second harmonic generation can simultaneously probe both the surface 
crystalline structure and the surface charge of Bi2Se3. Furthermore, we find that second harmonic 
generation using circularly polarized photons reveals the time-reversal symmetry properties of the 
system and is surprisingly robust against surface charging, which makes it a promising tool for 
spectroscopic studies of topological surfaces and buried interfaces. 



Electrons on the surface of a three-dimensional (3D) 
topological insulator [H-Q are predicted to exhibit exotic 
electrical properties including protection against local- 
ization from non-magnetic impurities P, [2] and photo- 
induced spin-polarized currents [1, Q • When interfaced 
with ordinary materials, these surface states are pre- 
dicted to evolve into new broken symmetry electronic 
phases with unconventional responses such as an anoma- 
lous half-integer quantum Hall effect 0, Q or topological 
superconductivity The recent discovery of 3D topo- 
logical insulators in Bii_2,Sba; 0, Bi2Se3 and related 
materials [9l-[Tl| has generated great interest to measure 
the symmetry and electrical properties of a single sur- 
face or buried interface. However transport techniques 
have difficulty separating the response from different sur- 
faces and the bulk [l^, and require deposition of 
contacts and gates that may perturb the electronic struc- 
ture. Optical techniques have been proposed as a contact 
free alternative that can be focused onto a single surface 
d, H, [H, [l^ or interface. However, experiments so far 
have been limited to the linear optical regime, which is 
again dominated by the bulk carrier response [la - llq . 

A contact free probe that is potentially surface or in- 
terface sensitive is second-order nonlinear optical spec- 
troscopy. In general, the electrical polarization of a ma- 
terial Pi [bj) has a dominant component linear in the driv- 
ing optical field Ej(u)) as well as weaker components pro- 
portional to higher powers of Ej{uj), where w is the op- 
tical frequency and the indices run through three spa- 
tial coordinates. Components that contain two powers 
of Ej{u}) are responsible for second harmonic genera- 
tion (SHG). For electric dipole processes, the polarization 
Pi = Xijk^jEk is obtained from a third rank susceptibil- 

ity tensor xljk that vanishes under inversion symmetry. 
Therefore dipole induced SHG is forbidden in bulk crys- 
tals with inversion symmetry and is only allowed at sur- 
faces or interfaces where inversion symmetry is necessar- 
ily broken [T9| . However, higher multipole bulk contribu- 
tions to SHG are still allowed even in inversion symmetric 
systems, which can obscure the surface dipole contribu- 
tion. Therefore it is unclear that SHG can be used as a 
surface or interface sensitive probe for a particular mate- 



rial. In this Letter, we generate SHG from Bi2Se3(lll) 
by utilizing the high peak photon flux of ultrafast lasers. 
By performing surface doping dependence measurements 
as a function of photon polarization and sample orienta- 
tion [Fig. 1(a)] we show that SHG is surface sensitive and 
is a simultaneous probe of both the surface crystalline 
structure and the surface Fermi level of Bi2Se3. Further- 
more, we find that circular dichroism SHG probes the 
time-reversal symmetry properties of the system and is 
robust against surface charging, which makes it a promis- 
ing tool for spectroscopic studies of topological surfaces 
and buried interfaces. 

Reflection SHG from the (111) surface of metallic 
Bi2Se3 bulk single crystals, which are bulk inversion sym- 
metric (space group D^^ fill), was measured using 80 fs 
laser pulses derived from a Ti: sapphire oscillator. The 
laser pulses have a center wavelength of 795 nm (1.56 
eV) and a repetition rate of 1.6 MHz, and are focused to 
a 20 /im spot on the sample at an incident angle of 45° 
with an incident power density of 0.63 kW/cm^, which 
is below the damage threshold. Specularly reflected pho- 
tons at the second harmonic energy (3.12 eV) were se- 
lectively isolated through spectral filtering and collected 
by calibrated photomultiplier tubes sensitive to 3.12 eV 
light. 

We first demonstrate that our measurements are con- 
sistent with SHG from the (111) surface where the crystal 
symmetry is reduced from Z)|^ to Cat,, which consists of 
an axis (z) of 3-fold rotational symmetry and three planes 
of mirror symmetry [Fig. 1(b)]. In general, the reflected 
surface SHG intensity I{2uj) is given by [20] 



I{2uj) ^Ax ]e.(2^)xge,(^)efc(u;)p/(c.) 



(1) 



where A is a constant determined by the experimental ge- 
ometry, e is the polarization of the incoming or outgoing 
photons, is the intensity of the incident beam and 

(2) 

Xijk ^ component tensor. Application of C^y sym- 

(2) 

metry reduces xljk '^'^^^ non-zero independent com- 
ponents Xxxx — Xxyy — Xyxyj Xzxx — Xzyyj Xxxz 

\ and Xzzz- By controlling the relative orientation be- 
tween the beam polarizations and crystal axes as shown 
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FIG. 1: (a) Schematic of the SHG experimental geometry. Surface and bulk regions are colored white and gray respectively, 
(b) Crystal structure of the Bi2Se3 (111) surface showing the topmost Se layer and underlying Bi layer, (c) Normalized SHG 
intensity I{2uj) from the (111) surface of Bi2Se3 measured as a function of azimuthal angle between the bisectrix (112) and 
the scattering plane. Measurements taken immediately after cleavage in Pin-Pout, (d) Sin-Pout, (e) Pin-Sout and (f) Sin-Sout 
incident and outgoing photon polarization geometries. Panels (g) to (j) shows analogous scans measured 200 minutes after 
sample cleavage. All data sets are normalized to the maximum intensity measured aX t = 200 mins in Pin-Pout geometry. The 
fewer data points in the early time scans is due to the need for faster data collection. Solid red lines are theoretical fits to 
Eqn.© [13]. 



in Fig. 1(a), components describing the in-plane and out- 
of-plane electrical responses can be isolated via: 

Ipp(2lj) = A X |a(3) - 0.025 a(i)cos(3(/))p 

/sp(2w) = ^ X |a(2) +0.016 a(i)cos(3(/))p 

Jps(2w) = A X 10.020 a(i)sin(30)p 

/S5(2a;) = ^ X 10.013 a(i)sin(30)|2 (2) 

where the first and second sub-indices of /(2a;) denote 
the input and output polarizations respectively, (j) is the 
angle between the scattering plane and the mirror (xz) 
plane of the crystal surface, a*-^-* = Xxxx is an in-plane 
response and a'^^ = 0.061 Xzxx and a^^^ = —0.007 Xxxz-^ 
0.096 Xzxx + 0.002 Xzzz are out-of-plane responses [21 1. 

Figure 1 shows the 4> dependence of SHG intensity 
from the (111) surface of Bi2Se3 measured under dif- 
ferent polarization geometries. The experiments were 
performed in air at room temperature and data were 
taken both immediately after and 200 minutes after the 
sample was cleaved. The time (t) dependence of the 
SHG patterns is related to surface charging as discussed 
later in the text. The SHG patterns exhibit a 6-fold 
or 3-fold rotational symmetry depending on whether S- 
or P-polarized output photons are selected respectively. 
Because ^-polarized light only has in-plane electric-field 
components, it is only sensitive to the in- plane response 
a*-^-' likely originating from anharmonic polarizability of 
Se-Se bonds [Fig. 1(b)], which have a 6- fold symmetric 
arrangement. On the other hand, because P-polarized 
light contains an electric-field component along z, it is 
sensitive to out-of-plane responses a^^' and a'^-' origi- 



nating from the Se-Bi bonds, which have a 3- fold sym- 
metric arrangement. A fit of equations ^ to the SHG 
patterns measured at a given time after cleavage using 
a single set of susceptibility tensor elements |2l| yields 
excellent agreement [Figs.l(c)-(f) & Figs.l(g)-(j)], which 
shows that the data is consistent with surface SHG from 
BizSeg. 

To prove that the surface contribution to SHG is domi- 
nant, we study its response to changes in the surface car- 
rier concentration. It is known from angle-resolved pho- 
toemission spectroscopy (ARPES) [11] that after cleav- 
age the surface electron concentration of Bi2Se3 increases 
monotonically over an hours long time scale before stabi- 
lizing. To study the effect of this surface electron doping 
on SHG we measure its time dependence immediately af- 
ter cleavage. Figure l(c)-(j) shows a clear change in the 
magnitude of all SHG patterns between t — and t = 200 
mins but no change in their rotational symmetry, which 
precludes a trigonal symmetry breaking atomic recon- 
struction. To investigate the possibility of surface atomic 
displacements that preserve overall trigonal symmetry 
such as a change of inter-layer distance, which would af- 
fect the in- and out-of-plane susceptibilities differently, 
we compare the complete time dependence of SHG in- 
tensities under different linear polarization geometries. 
Figure 2 shows that different combinations of suscepti- 
bility components all undergo the same monotonic in- 
crease by as much as 400% between t — and t ~ 50 
mins following cleavage, after which they saturate to a 
value that remains constant out to at least 600 mins, 
a trend highly consistent with the evolution of the sur- 
face Fermi level observed using ARPES ^2^] . To rule out 
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FIG. 2: Normalized SHG intensity along high symmetry di- 
rections of the (111) surface Brillouin zone of Bi2Se3 measured 
as a function of time after cleavage in air with (a) Pin-Pout, 
(b) Sin-Pout, (c) Pin-Sout and (d) Sin-Sout photon polariza- 
tion geometries. Insets in panels (c) and (d) show the time 
evolution of peak intensities measured starting 100 minutes 
after cleavage in air, prior to which the sample was not ex- 
posed to laser light. 



photo-induced effects, we performed time dependence ex- 
periments without exposing the sample to light for the 
first 100 mins after cleavage. In this case we find that the 
SHG intensities remain constant at the previous satura- 
tion values [insets Fig. 2(c) & (d)], which indicates that 
the initial rise in SHG intensity is purely a surface dop- 
ing effect and not a photo-induced effect [23l - [25j . Such 
large percentage changes in SHG intensity with surface 
electron doping suggest that SHG comes predominantly 
from electrons at the surface. 

Microscopically the increased surface carrier concen- 
tration affects SHG intensity through the perpendicular 
surface electric field £ that it creates. SHG is highly 
sensitive to a static electric field because it breaks in- 
version symmetry over the region that it penetrates, 
which is typically around 20 A for metallic Bi2Se3 (26j . 
Electric-field induced SHG is commonly observed in 
the context of metal/electrolyte interfaces, and is the- 
oretically described by a third order process Pi(2a;) = 

xf^ki^j{Q)ek{uj)ei{oj) 0,[2ll that acts in addition to x*^^- 

Because x^'^-' has the same symmetry constraints as x''^'' 
(2l| , the overall symmetry of the SHG patterns must re- 
main unchanged in the presence of £, with each tensor 
element in equation ([2|) simply being enhanced by the 
addition of a x*''^''^ tensor element. To demonstrate that 
SHG is sensitive to tuning the surface Fermi level we de- 
posit O2 on the surface, which has been shown to be an 
effective electron acceptor on Bi2Se3 (111) [IS]- Figure 
3 shows the time dependence of Ips{2lo) after cleaving 



(a) Surface molecular doping (b) SHG active regions 




FIG. 3: (a) Selected time dependent peak intensity curves 
measured in Pin-Sout polarization geometry after cleaving in 
O2. A similar decrease was also observed in other polariza- 
tion geometries, (b) Schematic of the energy evolution of the 
bulk conduction band minimum (Ec) and bulk valence band 
maximum (Ey) relative to the Fermi level (Ep) as a func- 
tion of distance to the O2 or (c) air covered surface. Typical 
Bi2Se3 samples show bulk conducting character {Ep > Ec) 
owing to Se vacancies 0, H^]. The gray, white and yellow re- 
gions correspond to the bulk, accumulation layer and surface 
of Bi2Se3. Schematic of the bulk (black) and surface (red) 
band structure is drawn near the surface region. 

in an O2 environment. We find that the SHG inten- 
sity exhibits an initial fast upturn reaching about 50% 
of the saturation intensity value in air, followed by a 
gradual decrease back to its initial value. These results 
show that electron transfer from the Bi2Se3 surface to 
O2 molecules takes place only after some finite £ has de- 
veloped, and that O2 can restore the surface back to an 
un-charged (£—0) state but cannot hole dope beyond this 
state. These behaviors are consistent with £ arising from 
the migration of negatively charged Se vacancies to the 
surface, which acts to lower the surface energy of the top- 
most Se layer. Assuming that charge transfer only takes 
place when an O2 molecule is adsorbed at a surface Se 
vacancy site, no additional hole doping will occur once all 
vacancies are occupied. Such a mechanism naturally ex- 
plains the slow time scale on which £ develops, as well as 
why the same behavior is observed regardless of whether 
the sample is cleaved in air or in ultra-high vacuum . 

Having established an optical method to probe the spa- 
tial symmetry and Fermi level of the surface electrons 
of Bi2Se3, we consider how SHG can be used to mon- 
itor time-reversal symmetry at a topological insulator 
surface. It has been proposed that new time-reversal 
symmetry broken phases can be measured through the 
differential index of refraction of right- (R) versus left- 
(L) circularly polarized light d, [l3,[ll|. However, it is 
known that the linear response of Bi2Se3 to circularly- 
polarized light is only sensitive to bulk electrons [l6| . In 
order to understand whether second harmonic versions of 
such experiments are feasible, it is necessary to measure 
the intrinsic second-order surface response of Bi2Se3 to 
circularly-polarized light. 

Second harmonic circular dichroism (CD), the differ- 
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FIG. 4: Normalized (/)-dependent SHG intensity patterns mea- 
sured 200 minutes after cleavage in air under (a) left-circular 
in (Lin) Sout and right-circular in (Rin) Sout, and (b) Lin- 
Pout and Rin-Pout photon polarization geometries. Solid lines 
are theoretical fits (see text), (c) Circular dichroism {Ir — Il) 
corresponding to data in panels (a) and (b). (d) Difference be- 
tween circular dichroism measured at 200 minutes after cleav- 
age and immediately after cleavage. 



ence in SHG yield using incident R- versus _L-polarized 
fundamental light, was measured from Bi2Se3 200 min- 
utes after cleavage in air. Figures 4(a) and (b) show 
that the SHG patterns obtained using R- and L-polarized 
incident light are clearly different for both P- and S- 
polarized output geometries, which are well described 
using equation ^ with the same set of susceptibility 
values fitted to Figs. 1(g) to (j). Although higher mul- 
tipole bulk contributions to SHG, which are allowed in 
inversion symmetric crystals, may contribute to CD due 
to an interference with the surface dipole radiation [l9| . 



such bulk terms are known to be greatly suppressed rel- 
ative to surface dipolar terms in the regime where the 
photon energy exceeds the bulk band gap [1^. Our ex- 
periment is in this regime and we have experimentally 
shown that the surface contribution is indeed dominant 
(Figs 1-3). While CD is generally non-zero, figure 4(c) 
shows that it vanishes when (j) is an integer multiple of 
60°, angles where the scattering plane coincides with a 
mirror plane of the (111) surface. Such zeroes are pro- 
tected by mirror symmetry because R- and L-polarized 
light transform into one another under mirror refiection 
about the scattering plane. Because a magnetization can 
break mirror symmetry, measuring departures from zero 
in CD along these specific values of (j) can be a sensitive 
probe of time-reversal symmetry breaking on the surface 
of Bi2Se3. Remarkably, we find that CD at a general 
(j) is insensitive to surface charging, as evidenced by the 
lack of measurable change as a function of time after 
cleavage [Fig. 4(d)]. This suggests that sensitive searches 
for time-reversal symmetry breaking induced CD may be 
carried out without the need for careful control of sur- 
face charging, which is an important and robust way of 
studying the physics of surface doped topological insu- 
lators or buried interfaces between topological insulators 
and ordinary materials [11, Q . 



The realization of surface second harmonic genera- 
tion from Bi2Se3 and the demonstration of its sensitivity 
to the surface Fermi level, crystal symmetry and time- 
reversal symmetry provides a novel contact-free probe of 
the properties of a single surface of a topological insula- 
tor. This will make it possible to probe the structural 
and electronic properties of a buried topological insula- 
tor interface, and to search for exotic broken symmetry 
interface phases. The use of ultrashort laser pulses will 
also make it possible to selectively study time-resolved 
non-equilibrium dynamics of a surface or interface fol- 
lowing an excitation, which for example can be applied 
to probe photo-induced spin-polarized currents. 
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